Energy and angular-momentum balance in wall-bounded superfluid turbulence 
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An example of steady-state wall-bounded quantum turbulence is provided by the turbulent vortex 
front. In a rotating cylinder filled with originally vortex-free superfluid 3 He-B, the front propagates 
axially and sets the superfluid in rotation. Our measurements of the front propagation at different 
angular velocities and temperatures down to 0.16T C show that besides the energy-dissipation mecha- 
nisms, the transfer of angular momentum is increasingly important for the turbulent dynamics in the 
T — > limit. We suggest a phenomenological model that describes the front velocity in terms of two 
effective mutual-friction parameters, which govern the energy and momentum balance, respectively. 
The residual values of these parameters in the T — !> limit differ by two orders of magnitude. In 
this limit, the efficient energy dissipation can be associated in traditional manner with the turbulent 
energy cascade toward small length scales, while the much less effective angular-momentum transfer 
causes decoupling of the superfluid from the rotating container. 

PACS numbers: 67.30hb, 47.15.ki, 47.20.-k, 67.25.dk 
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Introduction:-The new frontier of superfluid dynam- 
ics is the limit where the density of normal excitations 
tends to zero, on approaching zero temperature. There 
all traditional dissipation mechanisms in vortex flow 
rapidly vanish, extrapolations of the hydrodynamics from 
higher temperatures in terms of the superfluid Reynolds 
number predict increasing tendency towards turbulence, 
and the situation is reminiscent of what is known as the 
viscous anomaly in the turbulence of viscous fluids. This 
limit has been studied in both superfluid 4 He and 3 Hc-B 
with the focus on the dissipation rate, i.e., on the energy 
transfer from the superfluid to the normal component [lj- 
yj. Here we consider also the coupling of the superfluid 
to its container, i.e., the angular-momentum transfer to 
the superfluid component in a rapid increase of rotation 
in an axially symmetric container, the so-called spin-up. 

In classical viscous fluids, the spin-up response is dom- 
inated by surface interactions and characterized by the 
creation of the Ekman boundary layer Q. The details 
of the dynamics are affected by turbulent, nonlinear pro- 
cesses [5j. In the absence of quantized vortex lines, the 
superfluid component does not respond to the rotation of 
its container. This state in rotation is called the vortex- 
free Landau state. In the presence of vortices the transfer 
of angular momentum to the superfluid component is as- 
sisted by the mutual-friction force between the vortex 
lines and the normal- fluid excitations jf|. In superfluid 
4 He, the small vortex-core diameter of order 1 A empha- 
sizes the role of pinning of vortices at the container walls 
and vortex-boundary interactions 0, H[ . 

In 3 Hc-B contained in a smooth-walled cylinder, sur- 
face effects, such as vortex pinning, are strongly sup- 
pressed compared to superfluid 4 He, owing to the more 
than two orders of magnitude larger diameter of the vor- 
tex core. The spin-up responses at low temperatures, 
where the normal component consists of a dilute gas of 



ballistic thermal excitations, can be divided into three 
categories: First, in the absence of remanent vortices, the 
superfluid remains free of vortices and thus irrotational 
up to angular velocities of a few rad/s depending on the 
wall roughness and the cylinder radius, while the normal 
component follows the rotation of the container due to 
diffuse wall scattering. Almost no kinetic energy or an- 
gular momentum is transferred to the fluid due to the 
negligible density of the normal fluid. Second, if the con- 
tainer walls are initially covered with a sufficient amount 
of remanent vortices, the spin-up takes place as laminar 
solid-body-like motion determined by the mutual-friction 
coupling in the bulk liquid In this case, the trans- 
fer rate of both energy and angular momentum decrease 
exponentially with temperature, unless some residual ef- 
fects intervene, as observed in the laminar spin-down of 
Ref. 0. 

Here we study a third possibility, in which the spin-up 
starts from localized turbulent vortex formation. This 
localized turbulent burst, started by remnants pinned 
on a rough sintered heat exchanger surface, leads to the 
formation of a vortex front propagating axially along a 
smooth- walled long quartz cylinder 0] ■ The azimuthally 
precessing and axially propagating front motion corre- 
sponds to an axially inhomogencous spin-up of the su- 
perfluid component from rest (in the initial vortex-free 
Landau state at v s = 0) to a state with approximately 
solid-body rotation at a velocity v s = D, s z x r in the cylin- 
drical coordinates (r, <p, z). The front consists of a tur- 
bulent core with a size comparable to the cylinder radius 
and a bundle behind it which is approximately in solid- 
body rotation Across the core, in which the vortices are 
bent to the side wall, the angular velocity changes from 
zero to f2 s . This differential precession winds the vortex 
lines to a twisted configuration which extends from the 
front to the vortex bundle behind it 
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At sufficiently high temperatures, above 0.4 T c , mu- 
tual friction provides the coupling between the reference 
frame of the container and the rotating frame of the su- 
pcrfluid so that excluding the twist, the superfluid behind 
the front is at rest in the container frame. The coupling 
between the frames decreases with decreasing tempera- 
ture as the normal fluid density is reduced exponentially 
in the T —> limit. Bolometric measurements at 0.20T C 
supported by numerical simulations in Ref. [IH showed 
the density of the vortices behind the front to fall well 
below the equilibrium value and the front propagation to 
be followed by a slow laminar spin-up of the superfluid. 

Here we present new extensive measurements and anal- 
ysis of the decoupling in terms of the transfer rate of 
energy and angular momentum inferred from the prop- 
agation velocity of the turbulent vortex front both as a 
function of temperature and rotation velocity D, of the 
container. With adequate thermal contact between the 
superfluid sample and the heat exchanger, one can study 
the transfer properties down to the limit where the spin- 
up response becomes practically independent on temper- 
ature. 

Phenomenological model: At temperatures above 
0.5T C , front motion is laminar and its velocity is deter- 
mined by the mutual- friction dissipation a(T), 



Vf w a(T)ClR, 



(1) 



where f2 is the rotation velocity, and R the radius of the 
cylindrical container fl3| . At lower temperatures, where 
the dissipative mutual-friction parameter decreases as 

H 



*(T)=/?exp(-A/T), 



(2) 



the front motion starts to deviate from the laminar be- 
havior (JlJ due to the extra dissipation provided by tur- 
bulent processes [H . Here /3 is a pressure-dependent con- 
stant and A is the pressure- and temperature-dependent 
superfluid energy gap. In the simplest model of the front, 
the vortex state behind the thin front is in solid-body-like 
rotation with the angular velocity fl s w Q. 

At low temperatures, the helically twisted vortex bun- 
dle behind the front is in a q uasicquilibrium state at the 
reduced velocity Q s ^ f2 [12j. In analogy with the high- 
temperature behavior (Eq. [!}, with replaced by f2 s , 
one may suggest that 



V f « a eS (T)n s R, 



(3) 



where the effective mutual- friction parameter a e s(T) is 
the sum of a term proportional to the mutual-friction 
parameter and a turbulent-induced term, i.e., 



a e s{T) = Ca(T) + a tU rb- 



(4) 



Here C is a constant close to unity and for simplicity, 
aturb is assumed to be independent both of temperature 
and rotation velocity. 



Our decoupling model considers the angular- 
momentum balance of the rotating superfluid in 
terms of a course-grained hydrodynamical equation for 
the superfluid velocity v s , which is given by 



<9vs 

at 



V/x - v s x (V x v s ) = F Q 



(5) 



Here \x is the chemical potential. The two forces on the 
right-hand side are the mutual-friction force between the 
superfluid and the normal-fluid components with veloci- 
ties v s and v n , 

F„ = -a w x ((v s - v n ) x (V x v s )) , (6) 

and the force due to line tension (l5| . 

F A = -A(V x v») x (V x Cj) . (7) 

The unit vector Cj is directed along the vorticity V x v s . 
The line-tension parameter A depends on the intervortex 
distance I = y k/2CI and the vortex core diameter a as 

A = (k/47t) ln(£/a), (8) 

where £ ~ 0.1mm, a is comparable to the superfluid co- 
herence length ^ ps 70 nm at 0.5 bar liquid pressure, and 
K ~ 0.07mm 2 /s is the circulation quantum. 

In equilibrium, the global force balance dictates (F a ) = 
(Fx) for the azimuthal components F a and F\. With the 
simple estimate V x v s ~ 2Q S one obtains 



(F a ) ~ 2aO s (0 - n s )R. 



(9) 



The azimuthal line-tension force vanishes in the solid- 
body approximation. However, since the vortex state is 
twisted, the vorticity u is not perfectly parallel with z. 
The twist is proportional to Q s and has to be normalized 
by flR, i.e., V x Cj ~ fi s /(J!R) yielding 



(F X ) ^ 2\n 2 j(nR). 



(10) 



The solid-body-like angular velocity is obtained from the 
force balance equation 



aO s (Sl - fl s )R = \0%/(Q,R). 



(11) 



Due to the approximative nature of Eqs. (O and (p~0|) . we 
replace the mutual friction parameter a(T) in Eq. (11) 
with a new mutual- friction parameter a arn (T) for angular 
momentum, 



res • 



(12) 



where C am is a constant of the order of unity and the con- 
stant a res accounts for possible temperature-independent 
residual effects related to the angular-momentum trans- 
fer. Solving Eq. (fTTT) for Q s yields 



a am Cl + XR- 2 



n 2 



(13) 
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FIG. 1. The experimental setup. The He-B sample is con- 
tained in a 15-cm-long quartz cylinder with 6 mm inner di- 
ameter. The front motion is triggered by a rapid increase of 
the rotation velocity Q, which leads to the creation of vor- 
tices in the heat exchanger volume. The upward-propagating 
front is detected with two NMR pick-up coils, which are 9 cm 
apart. The quartz tuning forks at the bottom are used for 
thermometry. 

In Ref. [l2| this result was written in terms of the two 
Reynolds numbers Re Q f» l/a(T) and Rca ~ CtR 2 /X. 
In our case, the former controls the spin-up due to the 
coupling to the normal component, while the latter de- 
scribes the tendency of the superfluid to decrease the 
vortex number due to the vortex-line tension between 
the vortices in the front and the bundle behind it. Set- 
ting a am = a in Eq. (fT3")) , we get a simple interpolation 
formula in terms of the two Reynolds numbers |l2[ 



n s /n = (1 +Re Q /Rc A )~ 



(14) 



However here, making use of we fit the measure- 
ments to the more extended model 



Vf = (Cot + Oturb) 



Cam" 



(C am a + a [es )n + XR~ 2 



n 2 R. (15) 



Experimental techniques:-Our experimental setup 
(Fig.QJ consists of a long quartz cylinder filled with liquid 
3 He at 0.5 bar pressure. At its bottom end, the sample 
tube opens to a sintered-silver heat exchanger. Its rough 
surface ensures that vortices are formed there at low rota- 
tion velocity. In the experiment, we rapidly increase the 
rotation velocity from zero to the desired target velocity. 
This creates a turbulent burst at the bottom of the quartz 
tube and triggers an upward-propagating vortex front. 
Owing to the smooth walls, the critical velocity of vortex 
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FIG. 2. Measured front velocity Vf(T, Q) with logarithmic 
axes. The solid lines represent the velocity fitted to Eq. (|15l) . 
At high temperatures, the f2-dependence is linear, while in 
the zero-temperature limit the velocity tends to a quadratic 
dependence. The dashed lines correspond to Vf oc H and 
Vf oc tt 2 behavior for comparison. 



formation is as high as 1.8 rad/s in the quartz tube and 
thus at lower velocities, the volume above the front re- 
mains in the meta-stable vortex-free Landau state, which 
eventually is displaced by the front. 

The propagation is monitored with two NMR pick-up 
coils. At temperatures above 0.2T C , the arrival of the 
front can be detected by tracing the NMR signal at the 
so-called counterflow peak The arrival is seen as a 
rapid decrease of the NMR absorption. At the lowest ro- 
tation velocity of our measurements, fl 0.4 rad/s, this 
detection method loses sensitivity at around 0.23T C , es- 
pecially in the case of the coil at the top which operates 
at a smaller value of steady magnetic field. Even at our 
highest rotation velocity the counterflow peak becomes 
practically invisible below 0.18T C . There the front prop- 
agation can be monitored with the frequency shift of the 
magnon condensate in a magnetic trap in the middle of 
the pick-up coil, which depends strongly on the amount 
of local counterflow v n — v s [16]. The small axial pinch 
coils around the NMR pick-up coils in Fig. [T] are used to 
provide the localization of the trap in the axial direction. 
Temperature is measured from the resonance width of a 
quartz tuning fork oscillator, which depends on temper- 
ature as Af = 11460Hzexp(-1.776T c /T) in the ballistic 
regime of quasiparticle transport. 

After each measurement, the rotation is stopped and 
the sample is warmed up to T ~ 0.7T C for about one 



hour to allow the remanent vortices to annihilate 17 1 
before cooling down again for a new measurement. Al- 
ternatively, the sample can be warmed up above T c upon 
which it may be cooled down immediately. 

Results:-Figurc [2] shows the front velocity Vf(T, CI) 
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FIG. 3. The scaled front velocity Vf = Vf/(QR) as a func- 
tion of temperature at four different rotation velocities. The 
solid lines indicate the velocity fitted to Eq. (|15[) . In the 
zero-temperature limit the velocity tends to a temperature- 
independent but rotation-velocity-dependent constant value 
determined by the residual terms owb and a rcs . The inset 
shows the scaled velocity data at 29 bar from Ref. 0] for 
O ~ 1 rad/s. The solid line shows the predicted velocity from 
Eq. (|15[) using the same four parameter values obtained from 
the 0.5 bar measurements but with /3 and A replaced with 
their 29 bar values. 

as a function of rotation velocity in the range = 
0.4 — 1.7 rad/s at eight different temperatures ranging 
from 0.16 to 0.4T C . In this temperature interval the mu- 
tual friction parameter a changes by roughly three orders 
of magnitude. At 0.4T C the front velocity follows the 
single-vortex model of Eq. (fTJ) and the measured data 
points depend almost linearly on O. As expected from 
Eq. (fT5|) . the dependence gradually changes from lin- 
ear toward quadratic at temperatures below 0.3 T c . As a 
function of temperature, the scaled front velocity at four 
different rotation velocities is shown in Fig. [3] This plot 
demonstrates that Vf at constant O becomes temperature 
independent below 0.2T C as reported in Ref. 0- The new 
measurements reveal that the limiting scaled velocity in- 
creases with increasing O as expected from Eq. ([15)) . 

The fits in Figs. M and [3] use C 1 C am , c^turb and ct r cs 
as common fit ting parameters taking j3 from measure- 
ments in Ref. [lfjj. The obtained values are C « 0.52, 
C am w 1.33, a turb ~ 0.20, and a rcs « 0.0019. The fitted 
curves arc in reasonable agreement with the measured 
data and reproduce well the limiting behavior both in 
the high- and in the low-temperature limits. Addition- 
ally, small axial asymmetries in the precessing front and 
in the bundle behind it cause oscillations in the NMR 
signal, which allows us to measure O s directly, albeit 
in a limited range of temperatures slightly above 0.2 T c . 
These measurements qualitatively agree with Eq. (fT4)). 
Using the same four parameters in Eq. (|15|) produces a 
reasonable agreement with the data measured earlier at 
29 bar in a different quartz container. This is demon- 
strated in the inset of Fig. [3] Thus, it seems that the 




0.5 1 1.5 

n (rad/s) 



FIG. 4. Nonmonotonic features in the front velocity. The 
left panel shows a zoomed view of the scaled front velocity 
V{ = Vf/ttR at temperatures 0.28 and 0.30 T c illustrating the 
maxima at Q— 1-1.2 rad/s. On the right is a zoomed view of 
the temperature dependence for two different rotation veloci- 
ties showing plateaus at rotation-velocity-dependent temper- 
atures. The solid lines show the velocity from Eq. ()15[) with 
the fitting parameters given in the text. 

phenomcnological parameters have no strong dependence 
on the vortex-core diameter. 

The residual term of the angular-momentum transfer 
a res is about four times larger but still comparable to the 
corresponding residual term observed for laminar vortex 
motion in Ref. [13]. The larger value is expected, since it 
reflects the additional residual processes induced by the 
turbulence in the front. In the turbulent front motion, 
the vortex ends move along the sidewall of the cylinder, 
which supposedly enhances the surface interactions and 
increases angular-momentum transfer, when compared to 
the motion of predominantly parallel vortex lines sliding 
along, e.g., the flat top quartz end plate of the sample 
tube in laminar vortex flow flCj| . 

Importantly, a rcs is still a factor of hundred smaller 
than the residual term ctturb describing the energy trans- 
fer in the turbulent energy cascade. This means that 
the turbulent processes effectively enhance the energy 
transfer and consequently, the front velocity becomes 
practically temperature independent below 0.2 T c . The 
angular-momentum transfer, however, is not efficiently 
enhanced by turbulence. With vanishing mutual fric- 
tion, this leads to the observed quasiequilibrium solid- 
body-like rotation where the superfluid angular velocity 
O s < O depends strongly on rotation velocity. 

At the intermediate temperatures, the velocity behaves 
in a peculiar manner at high O, namely, Vf is nonmono- 
tonic with a temperature-dependent maximum at 0=1- 
1.2 rad/s as demonstrated in the left panel of Fig.|4] This 
is reflected as plateaus at rotation-velocity-dcpcndcnt 
temperatures ranging from 0.24-0.26T c for 1.4 rad/s to 
0.32-0.35T c for 0.5 rad/s as can be seen in the zoomed 
view of the temperature dependence in the right panel of 
Fig. |U A possible explanation is the bottleneck accumu- 
lation of the kinetic energy that leads to plateaus in the 
temperature range where the energy cascade approaches 
the intervortex distance I = yj k/ (20 s ) [2, 18. jl3|. This 
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assumption is supported by the fact that the value of a e g 
at the two plateaus decreases with increasing rotation ve- 
locity (see Fig. SI right). 



Discussion:-Our model qualitatively explains the ob- 
served front velocity in a wide region of angular veloci- 
ties (0.4-1.7 rad/s) and temperatures (0.16-0.4T C ). The 
origin of the observed residual term aturb has been at- 
tributed to the energy transfer in the turbulent energy 
cascade and to reflect the phenomenon of anomalous 
dissipation 0. Another recently suggested mechanism 
which could lead to temperature-independent dissipation 
in the zero-temperature limit is the leakage of quasipar- 
ticle excitations from overheated vortex cores due to the 
turbulence- induced Kelvin waves on vortex lines uOj . 



A new observation is the small value of the residual 
term a res in comparison to aturb- Consequently in the 
zero-temperature limit angular momentum transfer be- 
comes weak. Possibly surface interactions at the cylindri- 
cal wall, which have little effect at higher temperatures, 
play some role in the observed residual coupling. 



Contrary to earlier indications [21], these measure- 
ments show that turbulence is not the inevitable dynamic 
response in the T — > limit. New mechanisms help to 
promote laminar flow, like the partial decoupling from 
the applied rotation drive, followed by slow laminar re- 
covery to the equilibrium state. Secondly, as already 
known from measurements of the decay of vortex tan- 
gles in superfluid 4 He and in 3 He-B [22], vortex flow at 
T = is not dissipationless, but new sources of residual 
dissipation and friction are observed. The identification 
of these sources in both Bose and Fermi superfluids re- 
mains a task for the next near future. 
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